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ABSTRACT 

In the present work, we reconsider the idea of holographic dark energy. One of its key points is the 
formation of the black hole. And then, we propose the so-called "pilgrim dark energy" based on the 
speculation that the repulsive force contributed by the phantom-like dark energy (it) < —1) is strong 
enough to prevent the formation of the black hole. We also consider the cosmological constraints 
on pilgrim dark energy by using the latest observational data. Of course, one can instead regard 
pilgrim dark energy as a purely phenomenological model without any physical motivation. We also 
• briefly discuss this issue. 

O ! PACS numbers: 95.36.+X, 98.80.Es, 98.80.-k 

m 



60" 



CO 

> 

CN 
CO 

o 

O 

CN 



X 



* email address: haowei@bit.edu.cn 



2 



I. INTRODUCTION 



Since the discovery of the current accelerated expansion of our universe, dark energy has been one of 
the most active fields in physics and astronomy [l| . Of course, the simplest candidate of dark energy is 
a tiny positive cosmological constant, A. However, as is well known, it is plagued with the cosmological 
constant problem (the fine-tuning problem) and the cosmological coincidence problem. Therefore, many 
alternative dark energy models have been proposed. 

Recently, the so-called holographic dark energy (HDE) has been considered as an interesting candidate 
of dark energy, which has been studied extensively in the literature. For a quantum gravity system, the 
local quantum held cannot contain too many degrees of freedom, otherwise the black hole forms and then 
the quantum field theory breaks down. In the black hole thermodynamics 0, Q , there is a maximum 
entropy in a box of size L, namely the Bekenstein-Hawking entropy bound Sbh, which scales as the area 
of the box ~ L 2 , rather than the volume ~ L 3 . To avoid the breakdown of the local quantum field theory, 
Cohen et al. proposed the so-called energy bound, which is more restrictive than the entropy bound. 
If pa is the quantum zero-point energy density caused by a short distance cut-off, the total energy in a 
box of size L cannot exceed the mass of a black hole of the same size Q, namely p\L 3 <m 2 L, where 
rrip = (87rG) -1 / 2 is the reduced Planck mass. The largest IR cut-off L is the one saturating the inequality. 
Therefore, one has 

p A = Zc 2 m 2 L~ 2 , (1) 

where the numerical constant 3c 2 is introduced for convenience. In the literature, many IR cut-offs L 
have been considered, such as the Hubble horizon H~ x @, the particle horizon Rh = a J* dt/a — 
a Jq da/ (Ha 2 ) @,|43j], the future event horizon Rh = a dt/a — a da/ (Ha 2 ) Q, the Ricci scalar 
curvature radius, which is actually proportional to the causal connection scale of perturbations in the flat 
universe Rcc — (H +2H 2 )~ 1 / 2 (9[, the formal generalization of i?cc, namely (aH 2 + fiH)~ x I 2 (Tjj, the age 
of our universe T — da/ (Ha) [ll|, the conformal age of our universe r\ = f Q dt/a = J° da/(a 2 H) (l2l |. 
the radius of the cosmic null hypersurface [l3| , the so-called conformal- age-like length [l4 1 , where H = a/a 
is the Hubble parameter; a = (1 + z)^ 1 is the scale factor of our universe; z is the redshift; we have set 
ao = 1; the subscript "0" indicates the present value of the corresponding quantity; a dot denotes the 
derivative with respect to cosmic time t. 

In the present work, we reconsider the idea of holographic dark energy. One of its key points is 
the formation of the black hole. If we can prevent the formation of the black hole, the energy bound 
proposed by Cohen et al. ,4|] could be violated. If the repulsive force is strong enough, it might resist the 
matter collapse and then black hole does not form. So, what can contribute the strong repulsive force? 
Nowadays, it is well known that dark energy can contribute a repulsive force, since its equation-of-state 
parameter (EoS) w < —1/3. However, we will see below that the repulsive force contributed by the 
quintessence-like dark energy (w > — 1) is not strong enough to prevent the formation of the black hole. 
Therefore, we focus on the phantom-like dark energy (w < —1). Firstly, it is well known that everything 
will be completely torn up before our universe ends in the big rip caused by the phantom-like dark energy. 
Even the black hole will also be completely torn up. So, one can find that the repulsive force contributed 
by the phantom-like dark energy is strong enough to destroy the black hole. Secondly, in the famous 
work by Babichev et al. [l5[ , it is shown that accretion of phantom-like dark energy is accompanied with 
the gradual decrease of the black hole mass. Of course, this issue has not been completely settled by 
now. In Gao et al. claimed that the physical black hole mass may instead increase due to the 
accretion of phantom energy, but the Cosmic Censorship Conjecture will be violated. Later, Gonzalez 
and Guzman [l7l ] presented the full non-linear study of phantom scalar field accreted into a black hole, 
and found that the accretion of the phantom scalar field into the black hole can reduce its area down to 
50% within time scales of the order of few masses of the initial horizon. In [l8| , Sun also claimed that 
in the phantom dark energy universe the black hole mass becomes zero before the big rip is reached. 
On the other hand, Jamil and Qadir [l9| claimed that the phantom energy accretion contribute to a 
decrease of the mass of the primordial black hole. Recently, Sharif and Abbas |20| also found that mass 
of the black hole decreases due to phantom accretion. Thirdly, Harada et al. [21| claimed that there 
is no self-similar black hole solution in an universe with a stiff fluid or scalar field or quintessence. It 
is worth noting that Akhoury et al. [44j independently obtained the similar result from a very different 
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perspective. Chapline [22j claimed that black holes might not exist in our real world, because the negative 
pressure contributed by dark energy might prevent black holes from forming. Note that in [2l], [22|, [iij 
they only require dark energy violating the strong energy condition. Of course, they are wrong, because 
the repulsive force contributed by the quintessence-like dark energy (w > — 1) is not strong enough to 
prevent the formation of the black hole. This point can be seen in (23|, where by considering some 
specific McVittie solutions, Li and Wang showed that black holes can exist in a Friedmann-Robertson- 
Walker (FRW) universe dominated by dark energy. However, they require that the weak energy condition 
should be satisfied (although the other three energy conditions can be violated). Here, we would like to 
stress that the phantom-like dark energy does not satisfy the weak energy condition (in fact it violates 
all energy conditions) . Therefore, the conclusion of Li and Wang [23| cannot be applied to the phantom- 



like dark energy (w < —1). On the other hand, Rahaman et al. [24| claimed that they have found 
an exact solution of spherically symmetrical Einstein equations describing a black hole with a special 
type "phantom" energy source. Unfortunately, this special type "phantom" energy source used in [24j is 
characterized by p = —p in fact. Obviously, it is not the commonly called phantom energy (p < —p), and 
hence the conclusion of [24| cannot be applied here. 

Together with these three arguments mentioned above, it is reasonable to speculate that the repulsive 
force contributed by the phantom- like dark energy (w < —1) is strong enough to prevent the formation 
of the black hole. Of course, we admit that this speculation need further and solid proofs which are 
still absent. Anyway, we suggest considering what will happen to the argument of holographic dark 
energy if this speculation is true. In Sec.|Hl we propose the so-called "pilgrim dark energy" based on this 
speculation. Then, in Sec. IIII1 we consider the cosmological constraints on pilgrim dark energy by using 



the latest observational data. Some concluding remarks are given in Sec. IIV1 

II. PILGRIM DARK ENERGY 

Here, we propose the so-called pilgrim dark energy (PDE) based on the aforementioned speculation that 
the repulsive force contributed by the phantom-like dark energy (w < —1) is strong enough to prevent 
the formation of the black hole. If this speculation is true, the energy bound proposed by Cohen et al. Q 
could be violated, namely the total energy in a box of size L could exceed the mass of a black hole of the 
same size i.e., paL 3 >m 2 L . Therefore, the first property of pilgrim dark energy is 

PA>m 2 p L- 2 . (2) 

To implement Eq. ([2]), the simplest way is to set 

PA = Zn 2 m A - s L- s , (3) 

where n and s are both dimensionless constants. Thus, from Eqs. © and (O, p\ ~ m^ s L~ s >m 2 L~ 2 
leads to L 2 ~ s > nip" 2 = ( 2 ~ s , where £ p — m^ 1 — 1.616 x 10 -33 cm is the reduced Planck length, which is 
extremely short length in fact. Obviously, since L > £ p in general, it is required that 

s < 2 . (4) 

As mentioned above, the second requirement of pilgrim dark energy is to be phantom-like, namely 

iu A <-l. (5) 

Here, in order to obtain the EoS of pilgrim dark energy, we should choose a specific cut-off L. Of course, 
the simplest choice is the Hubble horizon L = H . In the epoch dominated by pilgrim dark energy, the 
Friedmann equation H 2 — > pa/ {3m 2 ) — n 2 m 2 ~ s H s leads to H —> const., which corresponds to wa — >• — 1. 
Therefore, our universe will end in a de Sitter phase, rather than big rip. In the matter-dominated epoch, 
H 2 - p m ~ a" 3 . So, p\ ~ H s ~ a" 3s / 2 . Thus, we have w A = -1 + s/2. From Eq. ©, s < is 
required. Let us consider the general case. Substituting pa oc H s into the energy conservation equation 
Pa + 3HpA (1 + wa) = 0, we have 



(6) 



4 



One might naively consider that H > in the late time, since wa < — 1. However, it is wrong. The total 
EoS Wtot = ^a^a — > — 1, where J7a is the fractional energy density of pilgrim dark energy. Therefore, 
H < instead holds in the late time. If H < is valid in the whole cosmic history, from Eqs. ([5]) and ©, 
it is required that 

s < . (7) 

So, there is no conflict between the matter-dominated epoch and the late time. w\ starts from —l + s/2 
in the matter-dominated epoch and goes asymptotically to —1 in the late time; u>a < — 1 always and 
it never crosses the phantom divide w = — 1 in the whole cosmic history. To be clearer, let us see a 
particular example of s = —2. In this case, substituting p m — p m od and Eq. ^ with L = H^ 1 into 
the Friedmann equation, we obtain 

H 2 = T-2 ("A + Pm) = n a r4H- a + |=f a~ 3 = c, + c 2 a' 3 . (8) 

Note that we consider a flat FRW universe which only contains pressureless matter and dark energy 
throughout this work. Requiring H 2 > 0, we can solve Eq. ([5} and obtain 
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H 2 = - c 2 a- 3 + JAc x + a- 6 . (9) 



When a — > 0, H 2 — > C2d~ 3 , which coincides with the one of matter-dominated epoch. When a — > oo, 
H 2 —> yfci = const., which coincides with the one of de Sitter phase in the late time. From Eq. ([5]), it 
is easy to see that H decreases when a increases; if is a monotonically decreasing function of a. Thus, 
H < indeed holds in the whole cosmic history. So, if s < 0, we have wa < —1 in the whole cosmic 
history. Note that if s = 0, from Eq. (O, it is easy to see that our pilgrim dark energy model reduces to 
the well-known ACDM model. Considering this point as well as Eqs. (Q} and ([7]), we allow 

s<0, (10) 

in order to satisfy the requirements © and ([5]), and include ACDM as a special case (although wa = — 1 
in this case). 

At first glance, one might consider that pilgrim dark energy model is a three-parameter model, which 
contains n, s and f2 m o as model parameters (where f2 m o is the present fractional energy density of 
pressureless matter). However, one of them is not independent in fact. Substituting p m — p m Qa~ 3 and 
Eq. ([3]) with L = H^ 1 into the Friedmann equation, namely 3m 2 H 2 = p\ + p m , we have 

„2 2-s tts-2 i PmO 3 -. ^ n. 

Introducing dimensionless E = H/Hq and h 2 = n 2 (m p / Hq) 2 ~ s , we recast Eq. (|1II) as 

h 2 E s ~ 2 + n m0 E- 2 a - 3 = 1 . (12) 

Requiring E{a = 1) = I by definition, we find that 

h 2 = i - n m0 • (13) 

So, n (or equivalently n) is not independent. Finally, the Friedmann equation becomes 

(i - a m0 ) e s - 2 + n m0 £- 2 (i + zf = \. (14) 

Obviously, there are only two free model parameters, namely f2 m o and s. From Eq. (|14p. one can obtain 
E(z) as a function of redshift z, if model parameters f2 m o and s are given. Note that E is a real number 
and E > is required by definition. 

Here, we would like to say some words before going further. As is shown above, pilgrim dark energy is 
based on the aforementioned speculation that the repulsive force contributed by the phantom-like dark 
energy (w < —1) is strong enough to prevent the formation of the black hole. In fact, this speculation 
is proposed just from some arguments, and it has no solid foundation so far. If one cannot agree the 
physical motivation of pilgrim dark energy presented here, we suggest that one can instead regard pilgrim 
dark energy as a purely phenomenological model without invoking any physical motivation. 
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FIG. 1: The 68.3% and 95.4% confidence level contours in the Q m o — s parameter space. The best-fit parameters 
are also indicated by a black solid point. This result is obtained by using the data of 580 Union2.1 SNIa alone. 



III. COSMOLOGICAL CONSTRAINTS ON PILGRIM DARK ENERGY 



In this section, we consider the cosmological constraints on pilgrim dark energy by using the latest 
observational data. 

At first, we use the observational data of Type la supernovae (SNIa) alone. Recently, the Supernova 
Cosmology Project (SCP) collaboration released the updated Union2.1 compilation which consists of 
580 SNIa [25|. The Union2.1 compilation is the largest published and spectroscopically confirmed SNIa 
sample to date. The data points of the 580 Union2.1 SNIa compiled in [25| are given in terms of the 
distance modulus Hobs{zi)- On the other hand, the theoretical distance modulus is defined as 

^th(zi) = 5log 10 D L (zi) + , (15) 
where /j,q = 42.38 — 51og 10 h and h is the Hubble constant Ho in units of 100 km/s/Mpc, whereas 

D L (z) = (l + z) f"-=^r, (16) 
Jo E(z;p) 

in which p denotes the model parameters. The x 2 from 580 Union2.1 SNIa is given by 

where a is the corresponding la error. The parameter /xq is a nuisance parameter, but it is independent 
of the data points. One can perform an uniform marginalization over /iq. However, there is an alternative 
way. Following [H, H3] , the minimization with respect to can be made by expanding the x 2 of Eq. (fT7|) 
with respect to n$ as 

xl(p) = A- 2tK>B + i4C, (18) 

where 

[fJ-obs(zi) ~ fJ>th(zi\ Mo = 0, p)] 2 
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B(p) = Z. ^-r^ ' c = <L 



Eq. (I18p has a minimum for /ig = B/C at 

^(p)=I(p)-^. (19) 

Since x 2 min = XL min obviously (up to a constant), we can instead minimize x 2 which is independent of 
/io- The best-fit model parameters are determined by minimizing the total \ 2 - When SNIa is used alone, 
we have \ 2 = x 2 which is given in Eq. (fl~9"|) . As in [2j],|28j], the 68.3% confidence level (C.L.) is determined 
by Ax 2 = X 2 — Xmm — l-0> 2-3 and 3.53 for n p = 1, 2 and 3, respectively, where n p is the number of 
free model parameters. Similarly, the 95.4% confidence level is determined by Ax 2 = X 2 ~ Xmm — 4-0, 
6.17 and 8.02 for n p = 1, 2 and 3, respectively. Here, we scan the 17 m o — s parameter space (note that 
as mentioned above, s < and < ri m o < 1 ar e required), and solve Eq. (fT4"| to obtain -E'(z) as a 
function of redshift z. Therefore, the corresponding \ 2 is on hand. Finally, we find that the best fit has 
Xmin = 562.226, and the corresponding best-fit parameters are ri m0 = 0.280 and s = —0.04. In Fig. [U 
we present the corresponding 68.3% and 95.4% C.L. contours in the O m o — s parameter space. It is easy 
to see that although the best-fit s is close to zero, there is a very big room for a significantly non-zero 
s in the 95.4% C.L. region. In fact, the viable s can extend to about —6.4 at 95.4% C.L., or —3.4 at 
68.3% C.L. The price to have a significantly non-zero s is a larger f2 m o- 

Next, we add the data from the observation of the large-scale structure (LSS). Here we use the distance 
parameter A of the measurement of the baryon acoustic oscillation (BAO) peak in the distribution of 
SDSS luminous red galaxies (29l. [30T|. which contains the main information of the observations of LSS. 
The distance parameter A is given by 



1 f Zb dz 1 2/3 



z b J E{z) 



(20) 



where z b = 0.35. In [30], the value of A has been determined to be 0.469 (n s /0.98)" - 35 ± 0.017. Here the 
scalar spectral index n s is taken to be 0.963, which has been updated from the WMAP 7-year (WMAP7) 




FIG. 2: The same as in Fig. [TJ except that this result is obtained by using the combined observational data of 
580 Union2.1 SNIa and the distance parameter A from BAO. 
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data 31]. Now, the total x 2 = xf, +X% where x; 2 is given in Eq. ([T9]). and x^i = (-4- .4 ;, s ) 2 /cr^. Again, 



we scan the f2 m o — s parameter space, and find that the best fit has X m in ~ 562.227, and the corresponding 
best-fit parameters are O m o = 0.278 and s — —0.008. In Fig. [21 we present the corresponding 68.3% and 
95.4% C.L. contours in the f2 m o — s parameter space. Comparing with Fig. [TJ it is easy to see that the 
contours are significantly shrunk. Although the best-fit s is very close to zero, there is still a room for 
a significantly non-zero s in the 95.4% C.L. region. In fact, the viable s can extend to about —1.05 at 
95.4% C.L., or -0.6 at 68.3% C.L. 
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FIG. 3: The same as in Fig. [TJ except that this result is obtained by using the combined observational data of 
580 Union2.1 SNIa, the distance parameter A from BAO, and the shift parameter R from WMAP7. 



Finally, we further add the data from the observation of the cosmic microwave background (CMB). 
Here we use the the shift parameter R, which contains the main information of the observations of the 
CMB [HH33. The shift parameter R of the CMB is defined by [H,[33[ 



R = a'/n / "7T7TT , (21) 



mO 



dz 



where the redshift of recombination z* = 1091.3 which has been updated in the WMAP7 data [31]. The 
shift parameter R relates the angular diameter distance to the last scattering surface, the comoving size 
of the sound horizon at and the angular scale of the first acoustic peak in CMB power spectrum of 
temperature fluctuations [H, [33|. The value of R has been updated to 1.725 ± 0.018 from the WMAP7 
data [3!]. Now, the total x 2 = X 2 + Xa + Xr: where Xr — (R ~ Robs) 2 /&%■ We scan the Q m0 - s 
parameter space, and find that the best fit has Xmin = 562.546, and the corresponding best-fit parameters 
are fi m o = 0.274 and s = 0.0. In Fig. [3J we present the corresponding 68.3% and 95.4% C.L. contours 
in the O m o — s parameter space. Comparing with Fig. [2J it is easy to sec that the contours are further 
shrunk. Although the best-fit s becomes zero, there is still a room for a non-zero s in the 95.4% C.L. 
region. In fact, the viable s can extend to about —0.56 at 95.4% C.L., or —0.29 at 68.3% C.L. 



IV. CONCLUDING REMARKS 



Some remarks are in order. Firstly, as mentioned above, pilgrim dark energy model reduces to the 
well-known ACDM model if s = 0. On the other hand, when we constrain pilgrim dark energy by 
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using the latest observational data, the best-fit parameter s is zero. However, this does not mean that 
pilgrim dark energy fails, because there is still a room for a non-zero s in the 95.4% C.L. region. For 
instance, pilgrim dark energy with s = — 1/2 is still viable. On the other hand, as is well known, ACDM 
model is plagued with the cosmological constant problem (the fine-tuning problem) and the cosmological 
coincidence problem. Therefore, pilgrim dark energy still deserves further investigations. 

Secondly, although wa < — 1 in the whole cosmic history, it will go asymptotically to —1 in the late 
time. So, our universe will end in a de Sitter phase, rather than big rip. This is an advantage of pilgrim 
dark energy model. 

Thirdly, in the present work, we have chosen IR cut-off L = H~ x for pilgrim dark energy, since the 
Hubble horizon H^ 1 is the simplest case. However, one can instead choose other IR cut-offs L for pilgrim 
dark energy, such as the particle horizon Rh = af Q dt/a = a da/ '(Ha 2 ) pfl, the future event horizon 
Rh = a J°° dt/a = a J™ da/(Hd 2 ) Q, the Ricci scalar curvature radius, which is actually proportional 
to the causal connection scale of perturbations in the flat universe Rcc = (H + 2H 2 )~ 1 / 2 Q, the formal 
generalization of Rcc, namely (aH 2 + fiH)~ x l 2 [l(|, the age of our universe T — da /(Ha) [ll|, the 

conformal age of our universe 77 = J^dt/a = Jq da/(a 2 H) [12], the radius of the cosmic null hypersur- 
face [l3|, the so-called conformal-age-like length This situation is very similar to holographic dark 
energy. It is possible to have completely new results (e.g. the cosmological constraints might be different) 
if L is changed. Of course, we welcome other authors to explore this possibility. 

Fourthly, as is shown above, pilgrim dark energy is based on the aforementioned speculation that the 
repulsive force contributed by the phantom- like dark energy (w < —1) is strong enough to prevent the 
formation of the black hole. In fact, this speculation is proposed just from some arguments, and it has 
no solid foundation so far. If one cannot agree the physical motivation of pilgrim dark energy presented 
here, we suggest that one can instead regard pilgrim dark energy as a purely phenomenological model 
without invoking any physical motivation. 

Finally, if we consider pilgrim dark energy as a purely phenomenological model without any physical 
motivation, the physical requirement s < can be given up. So, pilgrim dark energy can have a richer 
phenomenology. On the other hand, we also note that pilgrim dark energy can include other existing 
dark energy models as its special cases. For instance, if s = 1 and L = H^ 1 , pilgrim dark energy reduces 
to the so-called QCD ghost dark energy (34l-l37j. and Dvali-Gabadadze-Porrati (DGP) model [HI, [3!|. 
Of course, if s = 2 and choosing various L, pilgrim dark energy reduces to holographic dark energy |8|, 
(new) agegraphic dark energy [ll], [l2j], Ricci dark energy Q, and so on. If L = H^ 1 and s is free, pilgrim 
dark energy reduces to the modified DGP model (a dark energy) [i(| , and the modified holographic dark 
energy with IR infinite extra dimensions 41]. If L = H~ x , Rh, Rh, and s is free, pilgrim dark energy 
reduces to the so-called holographic cosmological "constant" derived from a generalized holographic 
relation between UV and IR cut-offs [13]. We anticipate that pilgrim dark energy with various L and 
a completely free parameter s (which is not necessary to be s < 0) is a rich mine. 
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